home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1995 April / macformat-023.iso / Shareware City / Developers / RLaB / rlib / rank.r < prev    next >
Encoding:
Text File  |  1994-09-20  |  626 b   |  28 lines  |  [TEXT/RLAB]
  1. //-------------------------------------------------------------------//
  2. //
  3. //  Syntax:    rank ( A )
  4. //        rank ( A , tol )
  5.  
  6. //  Description:
  7.  
  8. //  Compute the rank of the matrix A. Rank returns the number of
  9. //  singular values that are larger than:
  10. //        max( size(x) ) * norm(x,"2") * eps. 
  11.  
  12. //  If the user specifies tol, the the number of singular values
  13. //  larger than tol is returned.
  14.  
  15. //-------------------------------------------------------------------//
  16.  
  17. rank = function(x, tol)
  18. {
  19.   global (eps)
  20.  
  21.   s = svd(x);
  22.   if (!exist (tol))
  23.   { 
  24.     tol = max(size(x)) * norm(x,"2") * eps;
  25.   }
  26.   return sum(s.sigma > tol);
  27. };
  28.